Quantum response to time-dependent external fields 
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Abstract. Recently, explicit real time dynamics has been studied in various systems. These 
quantum mechanical dynamics could provide new recipes in information processing. We study 
quantum dynamics under time dependent external fields, and explore how to control the 
quantum state, and also how to bring the state into a target state. Here, we investigate a 
pure quantum mechanical dynamics, dynamics in quantum Monte Carlo simulation and also in 
quantum master equation. For the control magnetic states, operators which do not commute 
with magnetization are important. We study case of the transverse Ising model, in which 
we compare natures of thermal and quantum fluctuations. We also study the cases of the 
Dzyaloshinsky-Moriya interaction, where we find a peculiar energy level structure. Moreover we 
study the case of itinerant magnetic state, where we study the change from the Mott insulator 
to the Nagaoka ferromagnetic state. Effects of dissipation are also discussed. 



1. Introduction 

Recently, quantum dynamics has attracted interest in the field of information processing, i.e., 
data propagation by quantum cryptography and quantum computing pQ , and quantum annealing 
method [2]. In order to manipulate the quantum state, we need to understand how the quantum 
state changes under a change of external field. In general, when we study quantum dynamics, 
we describe a state of the system by a wavefunction which is generally a superposition of the 
eigenstates. Therefore, the quantum state covers all the configuration space (Hilbert space) and 
has an advantage in the data processing. In particular, the quantum annealing makes use of 
this advantage to find the desired state. Peculiar quantum properties have been found in the 
so-called single molecular magnets, such as Mn^p], Fes [I], Vi5[5], etc. These systems could 
be candidates for a storage of information. Here, we will study some properties of quantum 
fluctuation. 



2. Quantum fluctuation 

2.1. single spin case 

First, let us study the quantum fluctuation of a single spin system denned by 



H T w = -Ha z -Fa x , (1) 
where o~ a denotes the a-th component of the Pauli operator 

"--(Ji)- < 2 » 

Here, we take the z axis as the quantization axis. 

In the case V = 0, the magnetization of the z-component is a good quantum number and no 
quantum fluctuations exist. The eigenstates are given by 

«r*|+) = |+), and o~ z \ — ) = — | — ). (3) 

The eigenenergies are E\ = —H and E2 = H, which are depicted by the dotted lines in the 
Fig. Ha) and we call them "diabatic state" . When we set an initial state at a negative field (say 
the point denoted by the open circle in the Fig. Ha)) and change the field to a positive large 
field, the state follows the curve and simply comes to the point denoted by the open triangle. 
Here, the spin state does not change at all. 

If we add the term of V (the transverse field), the energy levels are given by the solid curve. 
In this case, the infinitesimally slow sweep of the field adiabatically leads the state to the point 
denoted by the closed circle. Here, the state is always in the ground state, which we call 
"adiabatic state". If we sweep the field with a finite speed, then some amount of the state 
is scattered to the excited state (the open triangle). This process is described by the Landau- 
Zener-Stiickelberg mechanism [HJ \7\, where the probability of staying in the ground state is given 

by 

where AE is the energy gap at the avoided level-crossing point {H = 0) and in the present 
case AE = 2T, and v is the velocity of the field dH(t)/dt, and AM is the difference of the 
magnetizations of the diabatic state (AM = 2). The magnetization processes M(t) for various 
sweep speeds are depicted in Fig. [D(b) 



2.2. Cooperative systems: The transverse Ising model 

The effects of quantum fluctuation on the ordering phenomena have been studied. The most 
typical model is given by the transverse Ising model [S] 

h = -jy j *M- r £°?> (5) 

<ij> i 

where (ij) denotes interacting pair sites. The transverse field causes tilt of spins in the classical 
case, and it causes spins to flip in the quantum case. Thus, this field reduces the correlation 
function of the z-components. In the case T = 0, there are two eigenstates | + +•••+) and 

I ) which are degenerate. We may consider these states as classical stable states. In 

the presence of T, the ground state is given by a linear combination of states. This state can be 
considered as a quantum mixing state of two classically stable states. Competition between the 
classical order due to the interaction JJ2<ij> ^t^j an d the quantum fluctuation due to TJ^i °~i 
causes a so-called quantum phase transition in the ground state. 
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Figure 1. (a) The eigenenergies as functions of the field(-ff). The dotted curves give those for 
r = 0, and the solid curves give eigenenergies for T = 0.1. (b) The time dependencies of the 
magnetization M{t) for various values of the sweeping speeds v = 0.005, 0.01, 0.02, 0.03, 0.05 and 
0.1. The bold solid curve denotes the adiabatic change. 

First we study the one-dimensional transverse Ising model. In Fig. EUa), we depict the 
eigenenergies as functions of H with V = 0.5 J {N = 6). In the case H = 0, a quantum phase 
transition takes place at J = T. Also at H = 0, the Hamiltonian is expressed in terms of fermion 
annihilation and creation operators, c q and cj as 



Using this dispersion relation, we can calculate the energies of low energy states. In Fig. EJ^b), 
we depict the size dependence of the energy gaps between the ground state and the excited 
states at T = 0.5J. We find that the energy gap between the ground state and the first excited 
state A.E01 becomes exponentially small which represents the symmetry breaking phenomenon 
below the critical point (T < J). On the other hand, the energy gap between the ground state 
and the second and third excited states, AEq2 and A£o3, respectively, stays almost constant. 
The time evolution of this magnetization under field sweep 



reflects the energy structure. Here, is the wavefunction in the dynamical process. When 
we slowly sweep the magnetic field H(t) the magnetization shows a stepwise magnetization 
process as a sequence of LZS transitions at the avoided level crossing points [9j. If we sweep 
fast, then the magnetization shows a size-independent shape, which we call "quantum spinodal 
decomposition" [10]. The process sweeping r at H = can be studied exactly [TT] . 

Although, in one dimension, the model has no ordered phase at finite temperature, in higher 
dimensions, the model has an order-disorder phase transition at a finite temperature for small 
values of r. At T = the system is the usual Ising model. Therefore, we have a phase diagram 
schematically depicted in Fig. [3l 

When we change parameters along the solid arrow (F = 0) quantum fluctuations are absent, 
and the thermal fluctuation induces the phase transition. On the other hand, when we change 




M(t) = <*(t)|£ !*(*)> 



(7) 




Figure 2. (a) Eigenenergies as functions of H of the one-dimensional transverse Ising model 
with r = 0.5J. (b) The energy gaps at H = and V = 0.5 J between the ground state and the 
first excited state (dots), the second excited state (circles) and the third excited state (crosses). 
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Figure 3. Schematic phase diagram of the transverse Ising model in high dimensions {D > 2). 



parameters along the dotted arrow, (T = 0), the thermal fluctuations are absent, and the 
quantum fluctuation induces the phase transition. 

3. Comparison of thermal and quantum fluctuations 

First, let us compare the natures of the thermal and quantum fluctuations. In order to describe 
both fluctuations, it is convenient to express the state by the density matrix 

M 

P = ^2 u >tj\ i )(j\- (8) 

i=l 

In the classical case, the diagonal elements represent the probabilities of the state i, i.e. 
P(i) = wa. In the quantum system, the off-diagonal elements give the quantum fluctuations. 
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Figure 4. Typical configurations of (a) thermal (T/T = 3/5) and (b) quantum fluctuations 
(T/r = 0.002) in the quantum Monte Carlo simulation. The red points denote up spins, and 
the blue ones down spins. 



The dynamics of the density matrix is given by the Bloch equation 

ih^- t p=[H,p\. (9) 

In order to find the density matrix for the equilibrium state, quantum Monte Carlo method 
can be used. Then, the density matrix is expressed by {d + l)-dimensional configuration using 
the Suzuki- Trotter decomposition which is a kind of path-integral representation of the density 
matrix [12]. The equilibrium density matrix of the model is expressed by 

p = e- l3H /Z, Z = Tre~ m . (10) 



In order to study the nature of fluctuation, we study matrix elements of the density matrix 
which correspond to snapshots of quantum Monte Carlo simulation. A d-dimensional transverse 
Ising model is mapped into a (d + l)-dimensional Ising model {erf} with anisotropic couplings 

j3J 1/3 

-f^real space = "T7> an< ^ -^imaginary time = ^ ln(tanh ^^ r -'^) - (H) 

If 0/M is small, i.e. at a high temperature, K Iea \ space is small and -/^imaginary time is large. Thus, 
the fluctuation occurs in the real space. On the other hand, at low temperatures, the fluctuation 
occurs in imaginary axis. In Fig. [H we depict typical configurations of (a) the thermal and (b) 
quantum fluctuations in the quantum Monte Carlo simulation. 



3.1. Reentrant phenomena 

In classical systems, we have pointed out that frustrated configurations cause non-monotonic 
developments of ordering. In some systems, the sign of the correlation function changes as a 
function of the temperature, i.e. an antiferromagnetic correlation appears at high temperatures 
while a ferromagnetic one at low temperatures. This kind of non-monotonic effect causes 
a reentrant phase transition where different types of phases appear successively when the 
temperature changes [13, 14]. There the distribution of density of states (entropy) takes an 




Figure 5. A typical example of reentrant type of the correlation function, (a) A frustrated 
lattice. The circles denote the spins o\ and 02, and the triangles denote the decoration spins 
Sfc. (b) Correlation function between the spin 1 and 2 as a function of the temperature 
(n = 4,8, 12, 16, and 20), and (c) correlation function between the spin 1 and 2 as a function of 
r. (n = 2,4,6 and 8). 



important role. Similar behavior has been observed in the transverse Ising model with the 
change of V [15J. There, the quantum fluctuation is affected by the frustration as well as 
the thermal fluctuation. A typical example of reentrant type V dependence of the correlation 
function is depicted in Fig. [5] for a frustrated lattice whose Hamiltonian is given by 

n 

^reentrant = J' <?l<?2 ~ J Y (°1 + ^Sfc, °i = ±1> and Si = ±1, (12) 

fc=l 

with J = 1 and J' = § J. 
3.2. Quantum annealing 

In order to find the ground state of a complicated system H.q, there have been proposed various 
methods. The so-called annealing method is one of the typical methods for this purpose. In 
the usual thermal annealing method, the thermal fluctuations produce possible candidates of 
states for the update. Monte Carlo method provides an ensemble of states for an equilibrium 
state by making use of a kind of Markov chain (master equation). If we use only local updates 
of the state, the realization of equilibrium ensemble is often difficult due to the frozen effects 
due to energy barriers and/or entropy barriers p3]. To avoid the freezing, there have been 
introduced various techniques such as multi-canonical Monte Carlo method [16J, temperature 
exchange Monte Carlo method [17] . etc. These methods employ a wide variety of fluctuations 
that cause the sampling to be more efficient, and accelerate to converge to the desired ensemble. 
The Swendsen-Wang algorithm [18] introduces a graphical representation to realize an efficient 
sampling in which cluster flips performed systematically. 

Quantum fluctuations have been also used to efficiently search for the ground state, a 
technique that is called quantum annealing [2] . There the transverse field generates the quantum 
fluctuations. In the limit V = oo, the state is the ferromagnetic state aligned to the x direction, 
which is the sum of all state 

\F X )= \<Ti,---,VN) (13) 

(Ti=±l,---,crjv=±l 



is the ground state. When we reduce T gradually, the ground state changes adiabatically to the 
ground state of the original system. Generally we believe that there is no level crossing during 
the process r — > 0. Thus when we gradually reduce T, the ground state moves to the ground 
state at r = 0. This is the idea of the quantum annealing and has been successfully applied 
to various systems. Recently, the proof of the convergence has been given |19| for the quantum 
case as well as the classical case |20| . 

Here let us consider mechanisms to find the ground state. In order to find ground state, the 
simplest method is the exact diagonalization. If the system has K states we need K 2 memory. 
In order to reduce the necessary memory several methods have been invented, e.g. the power 
method, Lanczos method, etc. Then we need the memory proportional to K. In the 1/2-spin 
system, K = 2 N . Thus, these methods correspond to a full search over all the states. On the 
other hand, the Monte Carlo method in the classical systems requires only memory of the order 
N. In the quantum Monte Carlo method, we need memory of the order N x N T , where N T is 
the number of points along the imaginary axis. Thus, both systems have advantage when N 
is large. As we saw in Fig. [3l beside the thermal annealing (the solid line) and the quantum 
annealing (the dotted line), there are many other paths to reach the ground state. It would be 
an interesting problem to find the optimal path starting from the point (T = T = oo). The 
quantum annealing method could be applied to very wide variety of systems. For example, 
applications to an classification into clusters [21] , and also to the variational Bayes inference [22] 
are being developed. 



4. Other types of quantum fluctuation 

The so-called single molecular magnets such as Mni2, Fes, V15, etc. have attracted interests 
because they show a sign of quantum dynamics through the discrete energy level structure. The 
degree of freedom inside the single molecules would be used as a storage of information. Above 
we have studied the quantum fluctuation in the transverse Ising model. In this section, we study 
other types of quantum fluctuations. In general, when an interaction 7i' does not commute with 
the order parameter M (the magnetization in the above case) 

[M,H']?0, (14) 

the magnetization is no longer good quantum number, and we say that the system has the 
quantum fluctuations. Then we find interesting properties in the dynamics of M under the 
change of parameters of the Hamiltonian. Here we introduce some of examples. 



4-1- Dzyaloshinsky-Moriya interaction in the triangle lattice 

As the important interaction which does not commute with the magnetization and causes a kind 
of quantum mixing is the Dzyaloshinsky-Moriya (DM) interaction 

H DM = E x Sj. (15) 

<ij> 

This interaction is characterized the vectors Dij and thus the vectors must be compatible with 
the symmetry of the lattice. Here, we study the effect of this interaction on a system consisting 
of three spins making an equilateral triangle. The Hamiltonian is given by 



H 3 = jY,Si- S i+ i -Y^Di-SiX S i+1 -Hj^Si, (16) 
where S4 = S\. Because of the symmetry, it is required that Df = Df = -Df and 
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Figure 6. Energy structure of an equilateral triangle lattice with DM interaction (D x = D z = 
0.2J) for (a) = 0°, (b) 9 = 45°and (c) = 90°, where H-D = cos6. 

where R is a matrix of the rotation of 120°. Here it should be noted that when H is parallel to D, 
the magnetization along the field commutes with the Hamiltonian and no adiabatic transition 
takes place (Fig. Ela)). On the other hand, if they are not parallel, avoided level crossing 
structures appear. In Fig. El we show the energy structures as a function of the field for the case 
the angle between H and D is 0°, 45° and 90°. Avoided level crossing structures also appear 
at the crossing of the states with M = 1/2 and M = 3/2. 

In Fig. WOb), we have the following interesting characteristic of the adiabatic change of the 
magnetization. We start from the ground state in a large positive H where the magnetization 
is almost +3/2. When we reduce the magnetization from it, the state follows the curve drown 
by a thick curve, and it goes to a level with M = —1/2 at large negative H. If we start from 
the ground state in a large negative H, it goes to the state of M = 1/2. This indicates that 
the adiabatic change does not follow the ground state. The double degeneracy of the 5=1/2 
states are characterized by the chirality, because the eigenstates of the states for the translation 
operation are e l2n ^ 3 and e l47r//3 . In the presence of the DM interaction the states with different 
chirality are not degenerate. The same property holds for other cases except the case of 6 = 0° 
(e.g. Fig. EJc)). Various interesting magnetization loops in a field cycling appear according to 
the structures of the adiabatic energy levels [23]. As we saw above, the DM interaction causes 
energy gaps which strongly depend on the direction of the field. In contrast, it is known that 
the hyperfine interaction between the electron spins and nuclear spins causes gaps which are 
independent of the field direction [25]. Recently, coherent quantum dynamics of driven Rabi 
oscillation was observed in V15 [26] . 

4-2. Particle conveyance by a trap potential 

When we consider the motion of particles, the operators of the momentum p and the position 
x do not commute. Thus, in the process of acceleration of the particle, there occurs various 
interesting quantum effects [27j . 

4-3. Realization of the Nagaoka-ferromagnetism by removal of an electron 
It is known that the total spin of the Hubbard model is zero in the half- filled bipartite lattices, 
while it takes the maximum value when an electron is removed. This mechanism is called 
"Nagaoka ferromagnetism" [28]. We can demonstrate an adiabatic change between these state 
if we add an extra lattice point and absorb an electron by a strong chemical potential in a 




Figure 7. (a) Magnetization process for v = 0.1,0.2, and v = 0.4. (b) Magnetization process 
for v = 0.001, v = 0.006, and v = 0.01. (From J. Phys. Soc. Jpn. 70 (2001) 3385.) 



magnetic field [29 



5. Dissipation effect 

Now, we discuss effects of environments on the adiabatic process, which we call 'magnetic Foehn 
effect' |30| . Here, we consider the effects in the simplest model (Eq. ([T])) to simulate the phonon- 
bottleneck effect found in a magnetic molecule V15 which belongs to this category [5]. 
We use a quantum master equation [311 [32] for various sweeping velocities. 

^ = -i [H, p it)] - A ([X, Rp (*)] + [X, Rp (t)]t) , (18) 

where X is a system operator through which the system couples with the bath, R is a bath 
operator, and A is the coupling constant between the system and the reservoir. From now 
on, we set parameters as T = 0.5, T = 1.0, and A = 0.001. In Fig. Eta), we present the 
magnetization curves for fast sweeping rates, v = 0.1,0.2, and 0.4. Here we clearly find that the 
magnetic plateau decreases when v increases, which is consistent with Eq. @. The increase of 
the magnetization after the plateau is a process of relaxation to the equilibrium state caused by 
the dissipation term. The dotted line there denotes the adiabatic curve of the magnetization. In 
the case of much slow sweeping rates, we again find the magnetic plateau as shown in Fig. EKh), 
although the LZS transition probability (Eq. ([4])) is almost one in these sweeping rates. Here, the 
sweeping rates are v = 2 x 10 -3 , 8 x 10~ 3 , and 2 x 10~ 2 . The dotted line denotes the isothermal 
curve of the magnetization at the present temperature. We should note that the magnetic 
plateau in this figure increases when v increases, which is opposite to the fast sweeping case. 

The quantum master equation brings the density matrix to that of the equilibrium of the 
system [32]. However, contact with the thermal bath causes modification of the system, which 
has been known as the Lamb shift. It has been shown that this effect is systematically taken 
into account by modifying the master equation [33]. 



6. Summary and Discussion 

Here we explored properties of quantum fluctuations and possible manipulations of them. 
Because the nature of quantum fluctuation is different from that of thermal fluctuation, it could 
play important roles in information processing. In particular, we expect that the nontrivial 
properties of quantum dynamics would provide new developments of information processing. 



Acknowledgments 

This work was partially supported by a Grant-in-Aid for Scientific Research on Priority Areas 
"Physics of new quantum phases in superclean materials" (Grant No. 17071011), and also by 
the Next Generation Super Computer Project, Nanoscience Program of MEXT. We also thank 
the supercomputer center, Institute for Solid State Physics, University of Tokyo for the use of 
the facilities. 



References 

[1] Nielsen M A and Chuang I L 2000 Quantum Computation and Quantum Information (Cambridge: Cambridge 
University Press) 

[2] Kadowaki T and Nishimori H 1998 Phys. Rev. E 58 5355; Das A and Chakrabarti B K 2005 Quantum 
Annealing And Related Optimization Methods (New York: Springer- Verlag) ; Santoro G E and Tosatti E 
2006 J. Phys. A: Math. Gen. 39 R393; Das A and Chakrabarti B K 2008 Rev. Mod. Phys. 80 1061 

[3] Thomas L, Lionti F, Ballou R, Sessoli R, Gatteschi D and Barbara B 1996 Nature 383 145; Friedman J R, 
Sarachik M P, Tejada J and Ziolo R 1996 Phys. Rev. Lett 76 3830 

[4] Sangregorio C Ohm T, Paulsen C, Sessoli R and Gatteschi D 1997 Phys. Rev. Lett. 78 4645; Wernsdorfer W 
and Sessoli R 1999 Science 284 233 

[5] Chiorescu I, Wernsdorfer W, Miiller A, Bogge H and Barbara B 2000 Phys. Phys. Lett. 84 3454 

[6] Zener C 1932 Proc. R. Soc. London Ser A137 697; Landau L 1932 Phys. Z. Sowjetunion 2 46; Stiickelberg 
E C G 1932 Helv. Phys. Acta 5 3207 

[7] Miyashita S 1995 J. Phys. Soc. Jpn. 64 3207; Miyashita S 1996 J. Phys. Soc. Jpn. 65 2734. 

[8] Chakrabati B K, Dutta A and Sen P 1996 Quantum Ising Phase Transitions in Transverse Ising Models 
(Heiderlberg: Springer- Verlag) 

[9] De Raedt H, Miyashita S, Saito K, Garcia-Pablos D and Garcia N 1997 Phys. Rev. B 56 11761 
[10] Miyashita S, De Raedt H and Barbara B unpublished 
[11] Dziarmaga J 2005 Phys. Rev. Lett. 95 245701 
[12] Suzuki M 1976 Prog. Theor. Phys. 56 1454 

[13] Miyashita S 1983 Prog. Theor. Phys. 69 714; Kitatani H, Miyashita S and Suzuki M 1986 J. Phys. Soc. Jpn. 

55 865; Miyashita S and Vincent E 2001 Eur. Phys. J B 22 203 
[14] Tanaka S and Miyashita S 2005 Prog. Theor. Phys. Suppl. 157 34 
[15] Tanaka S and Miyashita S unpublished 

[16] Berg B A and Neuhaus T 1991 Phys. Lett. B 267 249; Berg B A and Neuhaus T 1992 Phys. Rev. Lett. 68 9 
[17] Hukushima K and Nemoto K 1996 J. Phys. Soc. Jpn. 65 1604 

[18] Swendsen R H and Wang J S 1987 Phys. Rev. Lett. 58 86; Wolff U 1989 Phys. Rev. Lett. 62 361 
[19] Morita S and Nishimori H 2006 J. Phys. A: Math. Gen. 39 13903; Morita S and Nishimori H 2007 J. Phys. 
Soc. Jpn. 76 064002 

[20] Geman S and Geman D 1984 IEEE Trans. Pattern Anal. Mach. Intell. 6 721 
[21] Kurihara K, Tanaka S and Miyashita S unpublished 

[22] Sato I, Kurihara K, Tanaka S, Nakagawa H and Miyashita S unpublished 

[23] Miyashita S and Nagaosa N 2001 Prog. Theor. Phys. 106 533; De Raedt H, Miyashita S, Michielsen K and 

Machida M 2004 Phys. Rev. B 70 064401 
[24] Choi K Y, Matsuda Y H, Nojiri H, Kortz U, Hussain F, Stowe A C, Ramsey C and Dalai N S 2006 Phys. 

Rev. Lett. 96 107202 
[25] Miyashita S, De Raedt H and Michielsen K 2003 Prog. Theor. Phys. 110 889 
[26] Bertaina S, Gambarelli S, Mitra T, Tsukerblet B, Miiller A and Barbara B 2008 Nature 453 203 
[27] Miyashita S 2007 J. Phys. Soc. Jpn. 76 104003 
[28] Nagaoka Y 1966 Phys. Rev. 147 392 
[29] Miyashita S 2008 Prog. Theor. Phys. 120 785 
[30] Saito K and Miyashita S 2001 J. Phys. Soc. Jpn.70 3385 

[31] Weidlich W and Haake F 1965 Z. Phys. 185 30; Kubo R, Toda M and Hashitsume N 1985 Statistical Physics 
II (New York: Springer- Verlag) ; Louisell W H 1973 Quantum Statistical Properties of Radiation (New 
York: Wiley) 

[32] Saito K, Takesue S and Miyashita S 2000 Phys. Rev. E 61 2397 

[33] Mori T and Miyashita S 2008 Dynamics of the Density Matrix in Contact with a Thermal Bath and the 
Quantum Master Equation to appear in J. Phys. Soc. Jpn. 77 (Preprint arXiv:0810.0626 ) 



